Seismic Analysis and Design of Composite Shear Wall with Stiffened Steel Plate and Infilled Concrete

Several experiments are conducted to investigate the seismic behavior of composite shear walls because of their advantages compared to traditional reinforced concrete (RC) walls. However, the numerical studies are limited due to the complexities for the steel and concrete behaviors and their interaction. This paper presents a numerical study of composite shear walls with stiffened steel plates and infilled concrete (CWSC) using ABAQUS. The mechanical mechanisms of the web plate and concrete are studied. FE models are used to conduct parametric analysis to study the law of parameters on the seismic behaviour. The finite element (FE) model shows good agreement with the test results, including the hysteresis curves, failure phenomenon, ultimate strength, initial stiffness, and ductility. The web plate and concrete are the main components to resist lateral force. The web plate is found to contribute between 55% and 85% of the lateral force of wall. The corner of web plate mainly resists the vertical force, and the rest of web plate resists shear force. The concrete is separated into several columns by stiffened plates, each of which is independent and resisted vertical force. The wall thickness, steel ratio, and shear span ratio have the greatest influence on ultimate bearing capacity and elastic stiffness. The shear span ratio and axial compression ratio have the greatest influence on ductility. The test and analytical results are used to propose formulas to evaluate the ultimate strength capacity and stiffness of the composite shear wall under cyclic loading. The formulas could well predict the ultimate strength capacity reported in the literature.


Introduction
Steel-concrete composite shear walls are composed of web plates, stiffened plates, infilled concrete, and steel studs. The ultimate strength capacity and energy dissipation capacity of composite shear walls are higher than those of reinforced concrete (RC) walls. Researchers have been interested in the application of composite shear walls in buildings [1][2][3]. The relevant specifications for composite shear walls, such as seismic codes ASCE 7-10 [4] and AISC 341-10 [5], are formulated by permitting the use of composite steel plate shear wall (C-PSW) systems in earthquake zones.
Researchers have conducted experimental studies on composite shear walls without boundary walls. Nie et al. [6], Mydin [7], Wright [8], Wang [9], and Nie [10] showed that a composite shear wall has high ultimate strength capacity and good ductility. The failure mode was local buckling of the web plate and fracture failure of the corner of the wall. The design formula for the width-to-thickness ratio of the steel plate was proposed. Zhang et al. [11] and Zhang et al. [12] showed that more channels could weaken the ultimate strength capacity and stiffness of the wall but could enhance the ductility and energy dissipation capacity of the wall. Increasing the axial compression ratio had little effect on the yielding bearing capacity and ultimate strength capacity of the wall but increased the capacity of stiffness degradation. Finally, formulas for calculating the ultimate strength capacity and initial stiffness were proposed. However, the formulas could not include all the relevant parameters, and the mechanical mechanism was not explicit due to the limitation of the test. Thus, numerical analysis of the composite shear wall is necessary.
Researchers have studied FE models and proposed design formulas for composite shear walls. Nguyen et al. [13], Epackachi et al. [14], and Rafiei et al. [15] established FE models and verified their accuracy. The damage process of infilled concrete under a cyclic load and the proportion of the contribution of steel to the total shear force were analysed. Parameter analysis of the connector in the wall showed that more connectors could improve the bearing capacity of the steel plate and that changing the spacing of the connector could affect the failure mode of the steel plate. Wei et al. [16] studied the axial compression performance of composite shear walls. The influence of distance-to-thickness ratios on the failure mode was studied, and a formula for the axial compression capacity of a composite shear wall was proposed. The higher axial compression ratio of wall [17,18] is beneficial to restrain the internal concrete and improve the compressive strength of concrete, so the energy dissipation capacity of composite shear wall is enhanced. Increasing the thickness of the steel plate [19][20][21][22][23] can increase the stiffness and ultimate bearing capacity of the wall, as the hysteretic curve of the wall is plumper. Varma et al. [24] simulated shear walls with different aspect ratios. When the aspect ratio was between 0.6 and 3.0, the coupling effect of the moment and shear force was obvious. Specifications [25][26][27] define the formula for the shear capacity of composite shear walls. The formulas for the shear bearing capacity and the flexural bearing capacity were given, but the formula for the flexural-shear coupling was not provided.
In summary, the seismic performance of the CWSC is usually affected by factors such as the axial compression ratio and the shear span ratio. The scholars have carried out the relevant experimental studies on the composite shear wall, laying a foundation for the research on the seismic performance of the wall. The ultimate bearing capacity and lateral stiffness of composite shear walls are important parameters of seismic performance. However, most of the relevant researches are qualitative studies, and the influence of parameters on the ultimate bearing capacity and lateral stiffness are not quantitatively analysed. Although some studies have obtained the formulas for the ultimate bearing capacity and lateral stiffness, they are all based on the test results. The predictive effect with other cross-sectional forms are unknown. In addition, the formulas do not consider the influence factors of the number of channels, so it is not comprehensive enough.
In this paper, FE models are established with ABAQUS to simulate the seismic behaviour of composite shear walls under the cyclic loading and validated by tests. Figure 1 is a schematic diagram of composite shear wall. The mechanical mechanism, stress distribution and failure modes of the composite shear wall are researched by finite element analysis. A comprehensive parametric study is carried out to investigate the influence of parameters, including wall thickness, steel thickness, and shear span ratio, on the seismic behaviour. Formulas are proposed to predict the ultimate strength capacity and stiffness of composite shear walls and are validated by tests and parametric analysis. In addition, the other forms of composite shear walls are fitted by the formulas.

Sample Design
The author conducted experimental research on composite shear walls [12]. The steel plates were extended into the RC beam to prevent the samples from pulling out. The CWSC samples were constructed using double steel plates as the external component. The concrete was constructed as the infilled component. This experiment was based on the study of a shear wall in a super high-rise building [3], and five types of specimens were designed at a 1:5 scale. The CWSC samples were divided into several channels by the stiffened plates, as shown in Figure 2. The parameters of the samples are listed in Table 1. All the samples had a height of 1050 mm and a rectangular cross-section width of 700 mm. The thickness of the steel plate was 3 mm. The shear span ratio of all the samples was 0.75. The diameter of the stud was 6 mm, and the height of the stud was 20 mm.

Loading Programme and Test Setup
All the samples had the same axial force ratio (where n = 0.5), which was defined as Equation (1) [4].

Sample Design
The author conducted experimental research on composite shear walls [12]. The steel plates were extended into the RC beam to prevent the samples from pulling out. The CWSC samples were constructed using double steel plates as the external component. The concrete was constructed as the infilled component. This experiment was based on the study of a shear wall in a super high-rise building [3], and five types of specimens were designed at a 1:5 scale. The CWSC samples were divided into several channels by the stiffened plates, as shown in Figure 2. The parameters of the samples are listed in Table 1. All the samples had a height of 1050 mm and a rectangular cross-section width of 700 mm. The thickness of the steel plate was 3 mm. The shear span ratio of all the samples was 0.75. The diameter of the stud was 6 mm, and the height of the stud was 20 mm.

Sample Design
The author conducted experimental research on composite shear walls [12]. The steel plates were extended into the RC beam to prevent the samples from pulling out. The CWSC samples were constructed using double steel plates as the external component. The concrete was constructed as the infilled component. This experiment was based on the study of a shear wall in a super high-rise building [3], and five types of specimens were designed at a 1:5 scale. The CWSC samples were divided into several channels by the stiffened plates, as shown in Figure 2. The parameters of the samples are listed in Table 1. All the samples had a height of 1050 mm and a rectangular cross-section width of 700 mm. The thickness of the steel plate was 3 mm. The shear span ratio of all the samples was 0.75. The diameter of the stud was 6 mm, and the height of the stud was 20 mm.

Loading Programme and Test Setup
All the samples had the same axial force ratio (where n = 0.5), which was defined as Equation (1) [4].  where f c is the axial compressive strength of the concrete; f y is the yield strength of the steel plates; and A c and A s are the cross-sectional areas of the concrete and steel plate, respectively. Horizontal loading was controlled by the force and displacement [28]. In the forceloading phase, the horizontal forces were 300 kN and 600 kN, and loading was performed in one cycle. Horizontal loading was changed to displacement loading when the drift ratio reached 0.25%. The displacement increment was 2.625/0.25% (in mm/drift), and loading occurred in two cycles. The displacement increment was increased to 5.25/0.5% (in mm/drift) when the test displacement reached 21/2% (in mm/drift). The test was stopped either when the horizontal force dropped below 85% of the maximum strength or when a constant vertical axial force could not be maintained. The loading history is illustrated in Figure 3. where fc is the axial compressive strength of the concrete; fy is the yield strength of the steel plates; and Ac and As are the cross-sectional areas of the concrete and steel plate, respectively. Horizontal loading was controlled by the force and displacement [28]. In the forceloading phase, the horizontal forces were 300 kN and 600 kN, and loading was performed in one cycle. Horizontal loading was changed to displacement loading when the drift ratio reached 0.25%. The displacement increment was 2.625/0.25% (in mm/drift), and loading occurred in two cycles. The displacement increment was increased to 5.25/0.5% (in mm/drift) when the test displacement reached 21/2% (in mm/drift). The test was stopped either when the horizontal force dropped below 85% of the maximum strength or when a constant vertical axial force could not be maintained. The loading history is illustrated in Figure 3.   The test setup is shown in Figure 4. The specimen was laid between the top steel Lbeam and the bottom steel beam, and the RC beams were fixed to the bottom steel beam and the top steel L-beam. The bottom steel beam was anchored to the ground, and the top steel L-beam was connected by three actuators, one in the horizontal direction and the others in the vertical direction. On the top steel beam, four supports were arranged to prevent the out-of-plane deformation of the specimen during the test.  The test setup is shown in Figure 4. The specimen was laid between the top steel L-beam and the bottom steel beam, and the RC beams were fixed to the bottom steel beam and the top steel L-beam. The bottom steel beam was anchored to the ground, and the top steel L-beam was connected by three actuators, one in the horizontal direction and the others in the vertical direction. On the top steel beam, four supports were arranged to prevent the out-of-plane deformation of the specimen during the test.
where fc is the axial compressive strength of the concrete; fy is the yield strength of the steel plates; and Ac and As are the cross-sectional areas of the concrete and steel plate, respectively. Horizontal loading was controlled by the force and displacement [28]. In the forceloading phase, the horizontal forces were 300 kN and 600 kN, and loading was performed in one cycle. Horizontal loading was changed to displacement loading when the drift ratio reached 0.25%. The displacement increment was 2.625/0.25% (in mm/drift), and loading occurred in two cycles. The displacement increment was increased to 5.25/0.5% (in mm/drift) when the test displacement reached 21/2% (in mm/drift). The test was stopped either when the horizontal force dropped below 85% of the maximum strength or when a constant vertical axial force could not be maintained. The loading history is illustrated in Figure 3.   The test setup is shown in Figure 4. The specimen was laid between the top steel Lbeam and the bottom steel beam, and the RC beams were fixed to the bottom steel beam and the top steel L-beam. The bottom steel beam was anchored to the ground, and the top steel L-beam was connected by three actuators, one in the horizontal direction and the others in the vertical direction. On the top steel beam, four supports were arranged to prevent the out-of-plane deformation of the specimen during the test.  The FE model is composed of five parts, including the outer steel plate, steel studs, infilled concrete, concrete beams and reinforced cages. The outer steel plate is composed of a web plate and stiffened plate.
In order to simplify the calculation, the simple element type should be selected as far as possible. Since the thickness of the outer steel plate is less than 1/10 of its length and width, the stress change in the thickness direction can be ignored. So, the outer steel plate uses shell element (S4R). The change of stress in the thickness direction cannot be ignored in steel stud and concrete, because the sizes in three directions have little difference. The steel studs, infilled concrete and concrete beam use the solid element (C3D8R). The reinforced cage uses a beam element (T3D2), as shown in Figure 5. The assembled FE model is shown in Figure 5f.

Part and Element of the FE Model
The FE model is composed of five parts, including the outer steel plate, steel studs, infilled concrete, concrete beams and reinforced cages. The outer steel plate is composed of a web plate and stiffened plate.
In order to simplify the calculation, the simple element type should be selected as far as possible. Since the thickness of the outer steel plate is less than 1/10 of its length and width, the stress change in the thickness direction can be ignored. So, the outer steel plate uses shell element (S4R). The change of stress in the thickness direction cannot be ignored in steel stud and concrete, because the sizes in three directions have little difference. The steel studs, infilled concrete and concrete beam use the solid element (C3D8R). The reinforced cage uses a beam element (T3D2), as shown in Figure 5. The assembled FE model is shown in Figure 5f.

Contact of FE Model
In the FE model, the friction contact mode is used between the steel plate and concrete. For the normal contact behaviour, the "constraint enforcement method" in the default option is chosen, and "hard contact" in pressure-over closure is considered. For the tangential behaviour, the penalty in the friction formulation is chosen, and isotropic directionality is considered. The tangential friction coefficient is 0.6 [29], as shown in Figure  6a.

Contact of FE Model
In the FE model, the friction contact mode is used between the steel plate and concrete. For the normal contact behaviour, the "constraint enforcement method" in the default option is chosen, and "hard contact" in pressure-over closure is considered. For the tangential behaviour, the penalty in the friction formulation is chosen, and isotropic directionality is considered. The tangential friction coefficient is 0.6 [29], as shown in Figure 6a.
In the test, the steel studs were welded on the web plate. Thus, the steel studs are tied on the steel plate in the FE model, as shown in Figure 6b.
In the test, the reinforcement cage and the steel studs were poured with concrete. The reinforcement cage and the steel studs were fixed in the concrete. Thus, the reinforcement cage and the steel studs are embedded in the concrete beam and infilled concrete respectively, as shown in Figure 6c,d.  In the test, the steel studs were welded on the web plate. Thus, the steel studs are tied on the steel plate in the FE model, as shown in Figure 6b.
In the test, the reinforcement cage and the steel studs were poured with concrete. The reinforcement cage and the steel studs were fixed in the concrete. Thus, the reinforcement cage and the steel studs are embedded in the concrete beam and infilled concrete respectively, as shown in Figure 6c,d.

Boundary Conditions
In the test, the bottom beam was fixed to the ground through four bolts. The top beam is connected to the L-shaped beam through ten bolts, which limits the angle of the top beam and bottom beam. The boundary condition of the bottom beam is a fixed end constraint, and the boundary condition of the top beam is a sliding constraint, as shown in Figure 7. Therefore, six degrees of freedom are constrained at the bottom beam (i.e., U1 = U2 = U3 = UR1 = UR2 = UR3 = 0), and four degrees of freedom are constrained at the top beam (i.e., U3 = UR1 = UR2 = UR3 = 0).

Boundary Conditions
In the test, the bottom beam was fixed to the ground through four bolts. The top beam is connected to the L-shaped beam through ten bolts, which limits the angle of the top beam and bottom beam. The boundary condition of the bottom beam is a fixed end constraint, and the boundary condition of the top beam is a sliding constraint, as shown in Figure 7. Therefore, six degrees of freedom are constrained at the bottom beam (i.e., U1 = U2 = U3 = UR1 = UR2 = UR3 = 0), and four degrees of freedom are constrained at the top beam (i.e., U3 = UR1 = UR2 = UR3 = 0).   A reference point is set on the top of the model, which is coupled with the top surface. The axial compression force and lateral force are applied on the reference point. The loading process used in the model is the same as the tests, as shown in Figure 7.

Steel Constitutive Model
The constitutive model of steel is a fourfold linear isotropic strengthening model. The constitutive curve uses the uniaxial stress-strain curve of the test and converts to the true stress-strain curve [13]. The constitutive curve of the steel is shown in Figure 8. A reference point is set on the top of the model, which is coupled with the top surface. The axial compression force and lateral force are applied on the reference point. The loading process used in the model is the same as the tests, as shown in Figure 7.

Steel Constitutive Model
The constitutive model of steel is a fourfold linear isotropic strengthening model. The constitutive curve uses the uniaxial stress-strain curve of the test and converts to the true stress-strain curve [13]. The constitutive curve of the steel is shown in Figure 8. In the test process, the expansion of microcracks in the steel would cause a decrease in the ultimate strength capacity and stiffness of the material. Thus, the constitutive model of the steel must be considered the plastic damage.
The Wierzbicki [30] damage criterion considered that the excessive accumulation of plastic damage strain was the main cause of steel strength and stiffness degradation.
When the equivalent plastic strain pl  of the steel exceeds the equivalent plastic damage strain pl 0  , the internal steel is fractured and starts to damage. Wierzbicki [30] and Cook [31] proposed a steel plastic damage model and considered the equivalent plastic strain pl 0  to be associated with the degree of the triaxial stress  of steel. The degree of the triaxial stress  is defined as shown in Equations Equations (2)-(4). In the test process, the expansion of microcracks in the steel would cause a decrease in the ultimate strength capacity and stiffness of the material. Thus, the constitutive model of the steel must be considered the plastic damage.
The Wierzbicki [30] damage criterion considered that the excessive accumulation of plastic damage strain was the main cause of steel strength and stiffness degradation. When the equivalent plastic strain ε pl of the steel exceeds the equivalent plastic damage strain ε pl 0 , the internal steel is fractured and starts to damage. Wierzbicki [30] and Cook [31] proposed a steel plastic damage model and considered the equivalent plastic strain ε pl 0 to be associated with the degree of the triaxial stress η of steel. The degree of the triaxial stress η is defined as shown in Equations Equations (2)-(4).
where σ m and σ Mises represent the hydrostatic pressure and von Mises stress of steel, respectively. σ 1 , σ 2 , and σ 3 represent the principal stresses of steel in three directions. This paper adopts the simplified plastic damage model [32]. The equivalent plastic damage strain model of steel is defined in Equations Equations (5)- (7).
Materials 2022, 15, 182 where C 1 is the equivalent damage plastic strain of steel under pure shear (η = 0); C 2 is the equivalent damage plastic strain of steel under tension (η = η 0 ), which can be defined by the necking area A R ; n is the hardening coefficient of steel; and n = 0.226 in the FE model. η 0 = 1/3 [33]. The simplified curve ε pl 0 − η curve is shown in Figure 9.
where C1 is the equivalent damage plastic strain of steel under pure shear (η = 0); C2 is the equivalent damage plastic strain of steel under tension (η = η0), which can be defined by the necking area AR; n is the hardening coefficient of steel; and n = 0.226 in the FE model. η0 = 1/3 [33]. The simplified curve pl 0  − curve is shown in Figure 9. The value of steel damage is related to the plastic displacement when the steel achieves the damage state, which is represented by the damage factor Ds. The value of Ds ranges from 0 to 1. Ds = 0 indicates that no damage occurred to the steel, and the unloading stiffness don't degrade. Ds = 1 indicates that the steel is completely destroyed, and the unloading stiffness drops to 0, as shown in Figure 10. The value of steel damage is related to the plastic displacement when the steel achieves the damage state, which is represented by the damage factor D s . The value of D s ranges from 0 to 1. D s = 0 indicates that no damage occurred to the steel, and the unloading stiffness don't degrade. D s = 1 indicates that the steel is completely destroyed, and the unloading stiffness drops to 0, as shown in Figure 10. In this paper, the steel damage factor value D is adopted, as shown in Equation (8 where pl u and f u represent the plastic displacement and the ultimate displacement of the steel, respectively.

Concrete Constitutive Model
In the composite shear wall, the stress-strain curve of the infilled concrete is different from that of the plain concrete due to the constraint of the steel plate. The peak strain of the infilled concrete restrained by the steel tube is enhanced, and the descending stage is gentler because of the improved ductility [29]. The plastic damage model can be adopted to simulate the mechanical performance of concrete under cyclic loading in ABAQUS. The concrete constitutive model is calculated by the formulas presented in Concrete Infilled by Steel Tube Structure-Theory and Practice [35]. The formulas are shown in Equations (9) and (10). In this paper, the steel damage factor value D is adopted, as shown in Equation (8) [34]: where u pl and u f represent the plastic displacement and the ultimate displacement of the steel, respectively.

Concrete Constitutive Model
In the composite shear wall, the stress-strain curve of the infilled concrete is different from that of the plain concrete due to the constraint of the steel plate. The peak strain of the infilled concrete restrained by the steel tube is enhanced, and the descending stage is gentler because of the improved ductility [29]. The plastic damage model can be adopted to simulate the mechanical performance of concrete under cyclic loading in ABAQUS. The concrete constitutive model is calculated by the formulas presented in Concrete Infilled by Steel Tube Structure-Theory and Practice [35]. The formulas are shown in Equations (9) and (10).
The stress-strain curve of concrete under tensile stress is defined by Equation (11).
where x = ε ε p ; y = σ σ p , σ p = 0. 26(1.25 f c ) 2/3 ; and ε p = 4.31 × 10 −5 σ p . The Poisson's ratio of the concrete at the elastic stage is 0.2, and the elastic modulus is calculated by using Equation (12) [35]: Therefore, the constitutive curves of concrete under compression and tension are shown in Figures 11 and 12. As = Section area of the steel tube. Ac = Section area of the infilled concrete.
The Poisson's ratio of the concrete at the elastic stage is 0.2, and the elastic modulus is calculated by using Equation (12) [35]:   The main parameters of the plastic damage model of concrete in ABAQUS are as follows: the expansion angle Ψ = 38°, the flow potential function eccentricity ε = 0.1, the ratio of the initial equivalent biaxial compressive strength to the initial uniaxial compres- As = Section area of the steel tube. Ac = Section area of the infilled concrete.
The Poisson's ratio of the concrete at the elastic stage is 0.2, and the elastic modulus is calculated by using Equation (12) [35]:   The main parameters of the plastic damage model of concrete in ABAQUS are as follows: the expansion angle Ψ = 38°, the flow potential function eccentricity ε = 0.1, the ratio of the initial equivalent biaxial compressive strength to the initial uniaxial compres- The main parameters of the plastic damage model of concrete in ABAQUS are as follows: the expansion angle Ψ = 38 • , the flow potential function eccentricity ε = 0.1, the ratio of the initial equivalent biaxial compressive strength to the initial uniaxial compressive strength f b0 /f c0 = 1. 16, the ratio of the second stress invariant of the meridional deflection under tension and compression K c = 0.6, and the coefficient of viscosity µ = 0.001.
The concrete cannot heal completely after cracks appeared, and the unloading stiffness changes. The damage factor D c is used to describe the behaviour, as defined by Equation (13). As shown in Figure 13, the damage factor D c is defined based on the energy method [36]. Sc and S0 represent the shaded area and the area of triangle OAB, respectively. The damage factor of concrete is shown in Figures 11 and 12.

Meshing
Sensitivity analyses are carried out to determine the mesh size of steel plate and con-

Simplification of the Finite Element Model
The FE model with more elements requires a longer calculation time. Thus, the whole model can be converted to a half model to ensure the efficiency of calculation. The displacement in the Z direction and the rotation in the X and Y directions are constrained on the half surface. The half model and surface constraints are shown in Figure 15. The lateral force of the half model is reduced by half in the same loading displacement compared

Meshing
Sensitivity analyses are carried out to determine the mesh size of steel plate and concrete. The equivalent lengths of solid element ( 3 √ V element ) and shell element ( 2 √ V element ) are used to describe the mesh size. Figure 14 shows the strain history of local buckling at the bottom of the wall under different equivalent lengths. The strain values remain constant when the equivalent length range from 20 to 30 and 10 to 17.5, respectively. Considering calculation efficiency and accuracy, 3 √ V element = 30 and 2 √ V element = 17.5 are selected as the mesh size for finite element models. Sc and S0 represent the shaded area and the area of triangle OAB, respectively. The damage factor of concrete is shown in Figures 11 and 12.

Meshing
Sensitivity analyses are carried out to determine the mesh size of steel plate and con-

Simplification of the Finite Element Model
The FE model with more elements requires a longer calculation time. Thus, the whole model can be converted to a half model to ensure the efficiency of calculation. The displacement in the Z direction and the rotation in the X and Y directions are constrained on the half surface. The half model and surface constraints are shown in Figure 15. The lateral force of the half model is reduced by half in the same loading displacement compared with the full model. Therefore, the lateral force in the half model should be multiplied by 2 to compare with the test.

Simplification of the Finite Element Model
The FE model with more elements requires a longer calculation time. Thus, the whole model can be converted to a half model to ensure the efficiency of calculation. The displacement in the Z direction and the rotation in the X and Y directions are constrained on the half surface. The half model and surface constraints are shown in Figure 15. The lateral force of the half model is reduced by half in the same loading displacement compared with the full model. Therefore, the lateral force in the half model should be multiplied by 2 to compare with the test. As shown in Figure 16, the hysteretic curve and skeleton curve of the half model are similar to those of the overall model. Thus, the half model can simulate the hysteretic performance of the test. Half models are adopted to simulate the test in the following finite element analysis.

Force-Displacement Curve
To verify the accuracy of the FE model, the FE results of the hysteretic curves, skeleton curves and characteristic parameters are compared with those of the test. Figure 17 shows the comparisons of hysteretic curves and skeleton curves between the FE model and test. Figure 17 shows that the skeleton curves of FE model are similar to the test results at every loading level except CWSC-2. In addition, FE model accurately simulates the unloading path, especially at large horizontal displacements. Thus, the FE model can simulate the hysteretic curve of the composite shear wall. As shown in Figure 16, the hysteretic curve and skeleton curve of the half model are similar to those of the overall model. Thus, the half model can simulate the hysteretic performance of the test. Half models are adopted to simulate the test in the following finite element analysis. As shown in Figure 16, the hysteretic curve and skeleton curve of the half model are similar to those of the overall model. Thus, the half model can simulate the hysteretic performance of the test. Half models are adopted to simulate the test in the following finite element analysis.

Force-Displacement Curve
To verify the accuracy of the FE model, the FE results of the hysteretic curves, skeleton curves and characteristic parameters are compared with those of the test. Figure 17 shows the comparisons of hysteretic curves and skeleton curves between the FE model and test. Figure 17 shows that the skeleton curves of FE model are similar to the test results at every loading level except CWSC-2. In addition, FE model accurately simulates the unloading path, especially at large horizontal displacements. Thus, the FE model can simulate the hysteretic curve of the composite shear wall.

Force-Displacement Curve
To verify the accuracy of the FE model, the FE results of the hysteretic curves, skeleton curves and characteristic parameters are compared with those of the test. Figure 17 shows the comparisons of hysteretic curves and skeleton curves between the FE model and test. Figure 17 shows that the skeleton curves of FE model are similar to the test results at every loading level except CWSC-2. In addition, FE model accurately simulates the unloading path, especially at large horizontal displacements. Thus, the FE model can simulate the hysteretic curve of the composite shear wall.
To verify the accuracy of the FE model, the FE results of the hysteretic curves, skeleton curves and characteristic parameters are compared with those of the test. Figure 17 shows the comparisons of hysteretic curves and skeleton curves between the FE model and test. Figure 17 shows that the skeleton curves of FE model are similar to the test results at every loading level except CWSC-2. In addition, FE model accurately simulates the unloading path, especially at large horizontal displacements. Thus, the FE model can simulate the hysteretic curve of the composite shear wall. As shown in Figure 17b, the ultimate strength capacity of CWSC-2 in the test is lower than that of the FE model. The weld fractured when the lateral force increased to a certain value in the test. The steel plate could not resist as a whole, and the lateral force decreased.

Characteristic Parameters
The characteristic parameters include the initial stiffness, ultimate strength, and ductility. The Li method [37] and the Park method [38] were adopted to calculate the initial stiffness and yielding displacement, respectively, as shown in Figures 18 and 19. The ductility μ is defined as the ratio of the ultimate displacement to the yield displacement. As shown in Figure 17b, the ultimate strength capacity of CWSC-2 in the test is lower than that of the FE model. The weld fractured when the lateral force increased to a certain value in the test. The steel plate could not resist as a whole, and the lateral force decreased.

Characteristic Parameters
The characteristic parameters include the initial stiffness, ultimate strength, and ductility. The Li method [37] and the Park method [38] were adopted to calculate the initial stiffness and yielding displacement, respectively, as shown in Figures 18 and 19. The ductility µ is defined as the ratio of the ultimate displacement to the yield displacement.

Characteristic Parameters
The characteristic parameters include the initial stiffness, ultimate strength, and ductility. The Li method [37] and the Park method [38] were adopted to calculate the initial stiffness and yielding displacement, respectively, as shown in Figures 18 and 19. The ductility μ is defined as the ratio of the ultimate displacement to the yield displacement.  As shown in Figure 20, FE model mostly overestimates the ultimate strength within an error of 10%. The error in the ultimate strength of CWSC-2 exceeds 10%, because the welds fractured when the lateral force reached the ultimate strength [12]. The fractures are the reason why the test samples have a lower seismic resistance. The FE model could simulate the initial stiffness and ductility of the composite shear wall well within an error of 20%. The FE model underestimates the ductility of the composite shear wall. Based on the analysis above, the FE model could well simulate the seismic behaviour of composite shear walls.

Failure Damage
In the test loading process, the web plate experienced severe local buckling at different positions with increasing horizontal displacement. The FE model could simulate the macro phenomenon. Since the failure damages of the samples were similar, the failure phenomenon of CWSC-1 was taken as an example for comparison with the simulation.  As shown in Figure 20, FE model mostly overestimates the ultimate strength within an error of 10%. The error in the ultimate strength of CWSC-2 exceeds 10%, because the welds fractured when the lateral force reached the ultimate strength [12]. The fractures are the reason why the test samples have a lower seismic resistance. The FE model could simulate the initial stiffness and ductility of the composite shear wall well within an error of 20%. The FE model underestimates the ductility of the composite shear wall. Based on the analysis above, the FE model could well simulate the seismic behaviour of composite shear walls. As shown in Figure 20, FE model mostly overestimates the ultimate strength within an error of 10%. The error in the ultimate strength of CWSC-2 exceeds 10%, because the welds fractured when the lateral force reached the ultimate strength [12]. The fractures are the reason why the test samples have a lower seismic resistance. The FE model could simulate the initial stiffness and ductility of the composite shear wall well within an error of 20%. The FE model underestimates the ductility of the composite shear wall. Based on the analysis above, the FE model could well simulate the seismic behaviour of composite shear walls.

Failure Damage
In the test loading process, the web plate experienced severe local buckling at different positions with increasing horizontal displacement. The FE model could simulate the macro phenomenon. Since the failure damages of the samples were similar, the failure phenomenon of CWSC-1 was taken as an example for comparison with the simulation.

Failure Damage
In the test loading process, the web plate experienced severe local buckling at different positions with increasing horizontal displacement. The FE model could simulate the macro phenomenon. Since the failure damages of the samples were similar, the failure phenomenon of CWSC-1 was taken as an example for comparison with the simulation.

Distribution of the Lateral Force
The lateral force is transferred from the top to the bottom along the height of wall. The lateral force of the wall is constant along the height of the wall. However, the bending moment is variable along the height of the wall, the distribution of lateral force in the web plate and concrete are different. To research the distribution of lateral force at the different horizontal displacements, five positions of the lateral force in the web plate and concrete are extracted, as shown in Figure 22.    The lateral force is transferred from the top to the bottom along the height of wall. The lateral force of the wall is constant along the height of the wall. However, the bending moment is variable along the height of the wall, the distribution of lateral force in the web plate and concrete are different. To research the distribution of lateral force at the different horizontal displacements, five positions of the lateral force in the web plate and concrete are extracted, as shown in Figure 22.

Distribution of the Lateral Force
The lateral force is transferred from the top to the bottom along the height of wall. The lateral force of the wall is constant along the height of the wall. However, the bending moment is variable along the height of the wall, the distribution of lateral force in the web plate and concrete are different. To research the distribution of lateral force at the different horizontal displacements, five positions of the lateral force in the web plate and concrete are extracted, as shown in Figure 22.      When the drift ratio increases from 0.16% to 0.50%, the proportion of lateral force of the web plate decreases but increases in the concrete. The stresses in the steel and concrete increase rapidly in the elastic-plastic stage. The stress increases faster in the concrete than the steel. Then, with the increase of drift ratio, the proportion of lateral force of the web plate increases, but the concrete decreases. The stress of concrete has reached ultimate stress. Then, the capacity of the bearing lateral force in concrete decreases. The steel reaches the hardening stage. Then, the stress of steel increases with increasing drift ratio. The added lateral force is resisted by the web plate.    When the drift ratio increases from 0.16% to 0.50%, the proportion of lateral force of the web plate decreases but increases in the concrete. The stresses in the steel and concrete increase rapidly in the elastic-plastic stage. The stress increases faster in the concrete than the steel. Then, with the increase of drift ratio, the proportion of lateral force of the web plate increases, but the concrete decreases. The stress of concrete has reached ultimate stress. Then, the capacity of the bearing lateral force in concrete decreases. The steel reaches the hardening stage. Then, the stress of steel increases with increasing drift ratio. The added lateral force is resisted by the web plate. When the drift ratio increases from 0.16% to 0.50%, the proportion of lateral force of the web plate decreases but increases in the concrete. The stresses in the steel and concrete increase rapidly in the elastic-plastic stage. The stress increases faster in the concrete than the steel. Then, with the increase of drift ratio, the proportion of lateral force of the web plate increases, but the concrete decreases. The stress of concrete has reached ultimate stress. Then, the capacity of the bearing lateral force in concrete decreases. The steel reaches the hardening stage. Then, the stress of steel increases with increasing drift ratio. The added lateral force is resisted by the web plate.

Mechanical Mechanism of the Wall
The shear wall is composed of three parts, including the web plate, infilled concrete, and stiffened plate, as shown in Figure 25. In Section 4.4.1, the web plate and infilled concrete resist most of the lateral force. The stiffened plate has a low contribution to resist the lateral force. Thus, only the mechanism of the web plate and infilled concrete could be analysed.

Mechanical Mechanism of the Wall
The shear wall is composed of three parts, including the web plate, infilled concrete, and stiffened plate, as shown in Figure 25. In Section 4.4.1, the web plate and infilled concrete resist most of the lateral force. The stiffened plate has a low contribution to resist the lateral force. Thus, only the mechanism of the web plate and infilled concrete could be analysed.  The mechanical behaviours of the shear wall are generally similar, so CWSC-1 is taken as an example for study.
This section studies the mechanical mechanism of the components under the cyclic load through the maximum and minimum principal stresses. Figure 27

Mechanical Mechanism of the Wall
The shear wall is composed of three parts, including the web plate, infilled concrete, and stiffened plate, as shown in Figure 25. In Section 4.4.1, the web plate and infilled concrete resist most of the lateral force. The stiffened plate has a low contribution to resist the lateral force. Thus, only the mechanism of the web plate and infilled concrete could be analysed.  The mechanical behaviours of the shear wall are generally similar, so CWSC-1 is taken as an example for study.
This section studies the mechanical mechanism of the components under the cyclic load through the maximum and minimum principal stresses. Figure 27   The mechanical behaviours of the shear wall are generally similar, so CWSC-1 is taken as an example for study.
This section studies the mechanical mechanism of the components under the cyclic load through the maximum and minimum principal stresses. Figure 27 Figure 28a shows that the bottom section almost uniformly borne the compressive stress after axial force. With increasing lateral force, the normal stresses in the boundary parts of bottom section increase, but the middle part of bottom section is almost constant. Figure 28b shows that the shear stresses in the boundary parts of bottom section are almost equal to zero after the axial force is applied. With increasing lateral force, the normal stresses in the boundary parts of bottom section are almost constant, and the shear stress in the middle part of bottom section increases. The figures indicate that the mechanical mechanism in the bottom section of the wall is borne of the bending moment, shear force, and axial force. The bending moment is resisted by the boundary parts, and the shear force is resisted by the middle part.   Figure 28a shows that the bottom section almost uniformly borne the compressive stress after axial force. With increasing lateral force, the normal stresses in the boundary parts of bottom section increase, but the middle part of bottom section is almost constant. Figure 28b shows that the shear stresses in the boundary parts of bottom section are almost equal to zero after the axial force is applied. With increasing lateral force, the normal stresses in the boundary parts of bottom section are almost constant, and the shear stress in the middle part of bottom section increases. The figures indicate that the mechanical mechanism in the bottom section of the wall is borne of the bending moment, shear force, and axial force. The bending moment is resisted by the boundary parts, and the shear force is resisted by the middle part.   Figure 28a shows that the bottom section almost uniformly borne the compressive stress after axial force. With increasing lateral force, the normal stresses in the boundary parts of bottom section increase, but the middle part of bottom section is almost constant. Figure 28b shows that the shear stresses in the boundary parts of bottom section are almost equal to zero after the axial force is applied. With increasing lateral force, the normal stresses in the boundary parts of bottom section are almost constant, and the shear stress in the middle part of bottom section increases. The figures indicate that the mechanical mechanism in the bottom section of the wall is borne of the bending moment, shear force, and axial force. The bending moment is resisted by the boundary parts, and the shear force is resisted by the middle part.   Figure 28c shows that the middle section uniformly resist the compressive stress after the axial force is applied. With increasing lateral force, the normal stress in the middle section is almost constant. Figure 28d shows that with increasing lateral force, the shear stresses in the middle part of middle section increase rapidly, and the shear stresses in the boundary parts of middle section increase almost slowly. These figures indicate that the bending moment in the middle section of wall is equal to zero. The mechanical mechanism in middle section and middle part is resisting the shear force and the boundary part is resisting the axial force.
The shear stresses of the web plate mainly distribute around the stiffened plate and decrease at the stud. Because the studs connect the concrete and the web plate, part of the shear force transfer to the internal concrete through studs. Therefore, when the wall is subjected to lateral force, the shear force of the web plate between the studs is less than the stiffened plate. Another important reason is that the stiffened plate could prevent the web plate from buckling, so the web plate near the stiffened plate could carry more shear force, as shown in Figure 29.     Figure 28c shows that the middle section uniformly resist the compressive stress after the axial force is applied. With increasing lateral force, the normal stress in the middle section is almost constant. Figure 28d shows that with increasing lateral force, the shear stresses in the middle part of middle section increase rapidly, and the shear stresses in the boundary parts of middle section increase almost slowly. These figures indicate that the bending moment in the middle section of wall is equal to zero. The mechanical mechanism in middle section and middle part is resisting the shear force and the boundary part is resisting the axial force.
The shear stresses of the web plate mainly distribute around the stiffened plate and decrease at the stud. Because the studs connect the concrete and the web plate, part of the shear force transfer to the internal concrete through studs. Therefore, when the wall is subjected to lateral force, the shear force of the web plate between the studs is less than the stiffened plate. Another important reason is that the stiffened plate could prevent the web plate from buckling, so the web plate near the stiffened plate could carry more shear force, as shown in Figure 29.   Figure 28c shows that the middle section uniformly resist the compressive stress after the axial force is applied. With increasing lateral force, the normal stress in the middle section is almost constant. Figure 28d shows that with increasing lateral force, the shear stresses in the middle part of middle section increase rapidly, and the shear stresses in the boundary parts of middle section increase almost slowly. These figures indicate that the bending moment in the middle section of wall is equal to zero. The mechanical mechanism in middle section and middle part is resisting the shear force and the boundary part is resisting the axial force.
The shear stresses of the web plate mainly distribute around the stiffened plate and decrease at the stud. Because the studs connect the concrete and the web plate, part of the shear force transfer to the internal concrete through studs. Therefore, when the wall is subjected to lateral force, the shear force of the web plate between the studs is less than the stiffened plate. Another important reason is that the stiffened plate could prevent the web plate from buckling, so the web plate near the stiffened plate could carry more shear force, as shown in Figure 29.   As shown in Figure 31a, the minimum principal stress of concrete is almost vertical, which indicates that the concrete mainly resists the vertical force. As shown in Figure 31b,c, the normal stress of the middle section is almost unchanged with increasing drift ratio. Because the middle of wall mainly bears the shear As shown in Figure 31a, the minimum principal stress of concrete is almost vertical, which indicates that the concrete mainly resists the vertical force. As shown in Figure 31a, the minimum principal stress of concrete is almost vertical, which indicates that the concrete mainly resists the vertical force. As shown in Figure 31b,c, the normal stress of the middle section is almost unchanged with increasing drift ratio. Because the middle of wall mainly bears the shear As shown in Figure 31b,c, the normal stress of the middle section is almost unchanged with increasing drift ratio. Because the middle of wall mainly bears the shear force and axial force, the shear force of concrete is small. The vertical stress is almost constant.
The normal stress at the bottom section is divided into four sections, each of which increases from left to right. The concrete wall is separated into several concrete columns by the stiffened plates, each of which is independent and resists the lateral force alone. Thus, the vertical stress in the bottom section of each concrete column develops gradually from the left to right under the lateral force and axial force.

Parametric Analysis
To simulate the parameters of shear walls in actual high-rise buildings, parametric analysis is conducted based on the section size and design parameters of the shear walls in a proposed high-rise building [39]. The influence rules of the different parameters are studied, and formulas for the bearing capacity and stiffness of the composite shear wall are proposed. In order to study the hysteric behavior of model in parametric analysis, the cyclic load is applied to the wall. The model boundary condition is consistent with the test.
The parameters include the wall thickness, steel ratio, shear span ratio, axial compression ratio, concrete strength, steel strength and channel length-to-width ratio. A standard model is used as a reference, and then the parameters are changed to conduct parametric analysis. The standard model parameters are chosen as follows: the wall thickness is 1200 mm, the steel ratio is 5%, the shear span ratio is 0.5, the axial compression ratio is 0.5, the axial compressive strength of the concrete is 60 MPa, the yield strength of the steel is 345 MPa, and the length-to-width ratio of the channel is 2.
In the parametric analysis, the wall thickness ranges from 600 mm to 1500 mm, the steel ratio ranges from 3.3% to 8.3%, the shear span ratio ranges from 0.3 to 1.5, the axial compression ratio ranges from 0.2 to 0.6, the axial compressive strength of the concrete ranges from 40 MPa to 100 MPa, the yield strength of the steel ranges from 235 MPa to 420 MPa, and the length-to-width ratio of the channel ranges from 1 to 3.

Influence Rules of the Parameters
The hysteretic curve and skeleton curve are obtained by each model. The influence rules of key design parameters are studied by skeleton curve, including elastic stiffness, ultimate strength capacity, and ductility.

Influence on the Wall Thickness
As shown in Figure 32, the elastic stiffness and ultimate strength capacity of the wall enhance with increasing wall thickness, but improving the wall thickness has little effect on the ductility. Increasing the wall thickness could improve the section area of the wall. However, increasing the wall thickness only increases the thickness of the infilled concrete, which contributes slightly to the ductility of the composite shear wall.

Influence on the Steel Ratio
As shown in Figure 33, the elastic stiffness, ultimate strength capacity and ductility of the wall are enhanced with the increase of the steel ratio. Increasing the steel ratio could improve the flexural rigidity of the wall and the section area of the steel. The steel plate in the wall could resist most of the lateral force. The ultimate strength capacity and ductility are improved significantly with the increasing steel ratio.

Influence on the Steel Ratio
As shown in Figure 33, the elastic stiffness, ultimate strength capacity and ductility of the wall are enhanced with the increase of the steel ratio. Increasing the steel ratio could improve the flexural rigidity of the wall and the section area of the steel. The steel plate in the wall could resist most of the lateral force. The ultimate strength capacity and ductility are improved significantly with the increasing steel ratio.

Influence on the Steel Ratio
As shown in Figure 33, the elastic stiffness, ultimate strength capacity and ductility of the wall are enhanced with the increase of the steel ratio. Increasing the steel ratio could improve the flexural rigidity of the wall and the section area of the steel. The steel plate in the wall could resist most of the lateral force. The ultimate strength capacity and ductility are improved significantly with the increasing steel ratio.

Influence on the Shear Span Ratio
As shown in Figure 34, the elastic stiffness and ultimate strength capacity of the wall weaken but the ductility enhances with an increase in the shear span ratio. Increasing the shear span ratio changes the mechanical mode of the composite shear wall. In addition, the failure mode changes from relatively rigid shear failure to relatively flexible flexural failure. The elastic stiffness and bearing capacity enhance but the ductility weakens.

Influence on the Shear Span Ratio
As shown in Figure 34, the elastic stiffness and ultimate strength capacity of the wall weaken but the ductility enhances with an increase in the shear span ratio. Increasing the shear span ratio changes the mechanical mode of the composite shear wall. In addition, the failure mode changes from relatively rigid shear failure to relatively flexible flexural failure. The elastic stiffness and bearing capacity enhance but the ductility weakens.

Influence on the Steel Ratio
As shown in Figure 33, the elastic stiffness, ultimate strength capacity and ductility of the wall are enhanced with the increase of the steel ratio. Increasing the steel ratio could improve the flexural rigidity of the wall and the section area of the steel. The steel plate in the wall could resist most of the lateral force. The ultimate strength capacity and ductility are improved significantly with the increasing steel ratio.

Influence on the Shear Span Ratio
As shown in Figure 34, the elastic stiffness and ultimate strength capacity of the wall weaken but the ductility enhances with an increase in the shear span ratio. Increasing the shear span ratio changes the mechanical mode of the composite shear wall. In addition, the failure mode changes from relatively rigid shear failure to relatively flexible flexural failure. The elastic stiffness and bearing capacity enhance but the ductility weakens.  Compared to the stress distribution in different shear span ratios at the ultimate strength capacity (Figure 35), the Mises stress of web plate almost reaches ultimate strength in lower shear span ratio, and the web plate approximately reaches plane shear state. With the increase of shear span ratio, parts of the web plate cannot reach the ultimate strength, or even exceed the yield strength. The results indicate the web plate cannot completely develop plasticity when the higher shear span ratio wall reaches the ultimate strength capacity. Because the higher shear span ratio wall is controlled by the bending-shear failure mode, and the stress is small at the top of wall where the bending moment is small. Therefore, the shear span ratio of composite shear wall needs to be no more than 0.75. state. With the increase of shear span ratio, parts of the web plate cannot reach the ultimate strength, or even exceed the yield strength. The results indicate the web plate cannot completely develop plasticity when the higher shear span ratio wall reaches the ultimate strength capacity. Because the higher shear span ratio wall is controlled by the bendingshear failure mode, and the stress is small at the top of wall where the bending moment is small. Therefore, the shear span ratio of composite shear wall needs to be no more than 0.75.

Influence on the Axial Compression Ratio
As shown in Figure 36, the elastic stiffness and ultimate strength capacity of the wall enhance but the ductility weakens with an increase in the axial compression ratio. A larger axial compression ratio could increase the interaction between the steel and concrete, which improves the elastic stiffness and ultimate strength. However, the larger compression force tends to cause the shear wall damage rapidly after the ultimate strength and ductility weaken. The ductility weakens with an increase in the axial compression ratio.

Influence on the Axial Compression Ratio
As shown in Figure 36, the elastic stiffness and ultimate strength capacity of the wall enhance but the ductility weakens with an increase in the axial compression ratio. A larger axial compression ratio could increase the interaction between the steel and concrete, which improves the elastic stiffness and ultimate strength. However, the larger compression force tends to cause the shear wall damage rapidly after the ultimate strength and ductility weaken. The ductility weakens with an increase in the axial compression ratio.
strength, or even exceed the yield strength. The results indicate the web plate cannot completely develop plasticity when the higher shear span ratio wall reaches the ultimate strength capacity. Because the higher shear span ratio wall is controlled by the bendingshear failure mode, and the stress is small at the top of wall where the bending moment is small. Therefore, the shear span ratio of composite shear wall needs to be no more than 0.75.

Influence on the Axial Compression Ratio
As shown in Figure 36, the elastic stiffness and ultimate strength capacity of the wall enhance but the ductility weakens with an increase in the axial compression ratio. A larger axial compression ratio could increase the interaction between the steel and concrete, which improves the elastic stiffness and ultimate strength. However, the larger compression force tends to cause the shear wall damage rapidly after the ultimate strength and ductility weaken. The ductility weakens with an increase in the axial compression ratio.

Influence on the Axial Compressive Strength of Concrete
As shown in Figure 37, the elastic stiffness of the wall enhances with an increase in the axial compressive strength of the concrete. Increasing the axial compressive strength of concrete has little effect on the ductility and ultimate strength of the wall, and the elastic modulus improves. While increasing the ultimate strength capacity and ductility are relevant to increase the steel ratio. Thus, increasing the axial compressive strength of concrete has little effect on the elastic stiffness and ductility.
As shown in Figure 37, the elastic stiffness of the wall enhances with an increase in the axial compressive strength of the concrete. Increasing the axial compressive strength of concrete has little effect on the ductility and ultimate strength of the wall, and the elastic modulus improves. While increasing the ultimate strength capacity and ductility are relevant to increase the steel ratio. Thus, increasing the axial compressive strength of concrete has little effect on the elastic stiffness and ductility. As shown in Figure 38, the ultimate strength capacity of the wall enhances with increasing yield strength of steel. Increasing the yield strength of steel has little effect on the ductility and elastic stiffness. Increasing the yield strength of steel could not improve the elastic modulus and ductility of steel. Therefore, increasing the yield strength of steel could not improve the elastic stiffness and ductility of the wall. As shown in Figure 39, the elastic stiffness and ultimate strength capacity enhanced, but the ductility weakened with the increase in the length-to-width ratio of the channel. With increasing the channel length-to-width ratio, the mechanical mode of infilled concrete changed from multi-section columns to a whole wall. Therefore, increasing the length-to-width ratio of the channel could enhance the elastic stiffness and ultimate strength capacity of the wall but weaken the ductility.  Figure 38, the ultimate strength capacity of the wall enhances with increasing yield strength of steel. Increasing the yield strength of steel has little effect on the ductility and elastic stiffness. Increasing the yield strength of steel could not improve the elastic modulus and ductility of steel. Therefore, increasing the yield strength of steel could not improve the elastic stiffness and ductility of the wall.
As shown in Figure 37, the elastic stiffness of the wall enhances with an increase in the axial compressive strength of the concrete. Increasing the axial compressive strength of concrete has little effect on the ductility and ultimate strength of the wall, and the elastic modulus improves. While increasing the ultimate strength capacity and ductility are relevant to increase the steel ratio. Thus, increasing the axial compressive strength of concrete has little effect on the elastic stiffness and ductility. As shown in Figure 38, the ultimate strength capacity of the wall enhances with increasing yield strength of steel. Increasing the yield strength of steel has little effect on the ductility and elastic stiffness. Increasing the yield strength of steel could not improve the elastic modulus and ductility of steel. Therefore, increasing the yield strength of steel could not improve the elastic stiffness and ductility of the wall. As shown in Figure 39, the elastic stiffness and ultimate strength capacity enhanced, but the ductility weakened with the increase in the length-to-width ratio of the channel. With increasing the channel length-to-width ratio, the mechanical mode of infilled concrete changed from multi-section columns to a whole wall. Therefore, increasing the length-to-width ratio of the channel could enhance the elastic stiffness and ultimate strength capacity of the wall but weaken the ductility.  Figure 39, the elastic stiffness and ultimate strength capacity enhanced, but the ductility weakened with the increase in the length-to-width ratio of the channel. With increasing the channel length-to-width ratio, the mechanical mode of infilled concrete changed from multi-section columns to a whole wall. Therefore, increasing the length-towidth ratio of the channel could enhance the elastic stiffness and ultimate strength capacity of the wall but weaken the ductility.  Compared the stress distribution in different length-to-width ratio of the channel at the ultimate strength (Figure 40), the Mises stress is large in the middle part and small in the boundary part. The results indicate the wall is controlled by the bending-shear failure mode. With the decrease of length-to-width ratio of the channel, the area of shear stress area is larger, while the tension and compression stress areas are smaller due to the separation by stiffened plates. The stress distribution of web plate is more uniform when the length-to-width ratio of the channel decreases, and the performance of wall could be fully developed. Therefore, the length-to-width ratio of the channel in the composite shear wall needs to be no more than 1.2. Compared the stress distribution in different length-to-width ratio of the channel at the ultimate strength (Figure 40), the Mises stress is large in the middle part and small in the boundary part. The results indicate the wall is controlled by the bending-shear failure mode. With the decrease of length-to-width ratio of the channel, the area of shear stress area is larger, while the tension and compression stress areas are smaller due to the separation by stiffened plates. The stress distribution of web plate is more uniform when the length-to-width ratio of the channel decreases, and the performance of wall could be fully developed. Therefore, the length-to-width ratio of the channel in the composite shear wall needs to be no more than 1.2.

Formulation for Bearing Capacity and Stiffness Prediction
The composite shear wall is usually designed by considering two aspects: (1) the elastic stiffness K0 should be considered in elastic design, and (2) the yielding bearing capacity Fy, yielding section stiffness Ky, and ultimate strength capacity Fmax are studied in elasticplastic design.

Prediction of the Ultimate Strength
The following are assumed to identify the resistance of the composite shear wall: 1) The section of the wall remains planar after flexural bending. 2) The tension of the concrete

Formulation for Bearing Capacity and Stiffness Prediction
The composite shear wall is usually designed by considering two aspects: (1) the elastic stiffness K 0 should be considered in elastic design, and (2) the yielding bearing capacity F y , yielding section stiffness K y , and ultimate strength capacity F max are studied in elastic-plastic design.

Prediction of the Ultimate Strength
The following are assumed to identify the resistance of the composite shear wall: (1) The section of the wall remains planar after flexural bending. (2) The tension of the concrete on the tension side is ignored. The compressive strength of the concrete on the compression side is reached, and the yield strength of the steel plate is reached. (3) The restrained action of the steel plate on the concrete is considered. (4) The steel plate and concrete fully interact.
Based on the above assumptions, the section stress distribution of the shear wall under horizontal and axial compression forces is shown in Figure 41.
on the tension side is ignored. The compressive strength of the concrete on the compression side is reached, and the yield strength of the steel plate is reached. 3) The restrained action of the steel plate on the concrete is considered. 4) The steel plate and concrete fully interact.
Based on the above assumptions, the section stress distribution of the shear wall under horizontal and axial compression forces is shown in Figure 41.  According to the resultant force equilibrium equation, Equations (14) and (15) could be determined. M f cc = f c · (1 + 1.8ξ sc ) f cc = Axial compressive strength of the constrained concrete [40]. l c = Distance from the constrained concrete to the neutral axis. l yc = Distance from the compression steel plate to the neutral axis. l yt = Distance from the tension steel plate to the neutral axis.
As the composite shear wall is subject to the coupling action of axial force, shear force and bending moment, the interaction of three factors needs to be considered when calculating the ultimate strength [41]. The formula used to calculate the ultimate strength capacity of the composite shear wall is proposed in Equation (16).
N = Axial force applied at the shear wall. M = Bending moment applied at the shear wall. V = Horizontal shear force applied at the shear wall. N u = Axial bearing capacity of the shear wall, determined by Equation (17). M u = Flexural bearing capacity of the shear wall, determined by Equation (15). V u = Ultimate strength capacity of the shear wall, determined by Equation (18) Figures 42 and 43. The formula could fit the test with an error less than 10%, and fit the parametric analysis with an error less than 20%. Thus, the formula could accurately predict the ultimate strength capacity.     Figures 42 and 43. The formula could fit the test with an error less than 10%, and fit the parametric analysis with an error less than 20%. Thus, the formula could accurately predict the ultimate strength capacity.   The specification formulas used to calculate the ultimate bearing capacity, including the JGJ/T 380-2015 [27], JGJ 3-2010 [43], AISC 341-16 [44], and AISC N690-12 [25], compared with the test and parametric analysis, as shown in Figure 44a,c, show that the values in the test and parametric analysis are higher than the values obtained by the formulas and exceed 20% error, indicating that the ultimate bearing capacity is seriously underestimated by the JGJ/T 380-2015, AISC 341-16, and AISC N690-12. Figure 44b shows that the values in the test are lower than values obtained by the formula, within the 20% error. Most values obtained in the parametric analysis are higher than the those obtained with the formulas and exceed 20% error, indicating that the ultimate bearing capacity obtained in the parametric analysis is underestimated by JGJ 3-2010. Therefore, the formula in this paper could predict the ultimate bearing capacity of CWSC more accurately. timated by the JGJ/T 380-2015, AISC 341-16, and AISC N690-12. Figure 44b shows that the values in the test are lower than values obtained by the formula, within the 20% error. Most values obtained in the parametric analysis are higher than the those obtained with the formulas and exceed 20% error, indicating that the ultimate bearing capacity obtained in the parametric analysis is underestimated by JGJ 3-2010. Therefore, the formula in this paper could predict the ultimate bearing capacity of CWSC more accurately.

Prediction of the Yielding Bearing Capacity
The same method is used to calculate the yielding bearing capacity of the composite shear wall. The formula used to calculate the yielding bearing capacity of the composite shear wall is proposed in Equation (19).
The symbols in Equation (19) Figures 45 and 46. The formula could fit the test with an error less than 20%. The values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. Thus, the formula could accurately predict the yielding bearing capacity.

Prediction of the Elastic Stiffness
According to previous research [12], the elastic stiffness K0 of the composite shear wall is composed of the flexural stiffness Km and shear stiffness Kv. The formula for elastic stiffness is shown in Equations (20)- (22). Equation (19) is used to fit the test and FE model results. The test and parametric analysis results are compared with the formula results, as shown in Figures 45 and 46. The formula could fit the test with an error less than 20%. The values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. Thus, the formula could accurately predict the yielding bearing capacity.

Prediction of the Elastic Stiffness
According to previous research [12], the elastic stiffness K0 of the composite shear wall is composed of the flexural stiffness Km and shear stiffness Kv. The formula for elastic stiffness is shown in Equations (20)- (22).

Prediction of the Elastic Stiffness
According to previous research [12], the elastic stiffness K 0 of the composite shear wall is composed of the flexural stiffness K m and shear stiffness K v . The formula for elastic stiffness is shown in Equations (20)- (22). G sc A sc = Shear stiffness of the shear wall, where G sc A sc = q e (G s A s + G c A c ). q e represents the shear stiffness reduction coefficient, as proposed by Equation (23).
where ξ represents the length-to-width ratio of the channel, n represents the axial compression ratio, and λ represents the shear span ratio.
A comparison of the test and parametric analysis results with the formula results, as shown in Figures 47 and 48. The formula could fit the test with an error less than 10%. The values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. Thus, the formula could accurately predict the elastic stiffness of wall. pression ratio, and λ represents the shear span ratio.
A comparison of the test and parametric analysis results with the formula results, as shown in Figures 47 and 48. The formula could fit the test with an error less than 10%. The values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. Thus, the formula could accurately predict the elastic stiffness of wall.

Prediction of the Secant Stiffness of the Yield Point
The same method is used to calculate the secant stiffness of the yield point. The formula is proposed in Equations (24)- (27).  where ξ represents the length-to-width ratio of the channel, n represents the axial compression ratio, and λ represents the shear span ratio. A comparison of the test and parametric analysis results with the formula results, as shown in Figures 47 and 48. The formula could fit the test with an error less than 10%. The values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. Thus, the formula could accurately predict the elastic stiffness of wall.

Prediction of the Secant Stiffness of the Yield Point
The same method is used to calculate the secant stiffness of the yield point. The formula is proposed in Equations (24)- (27).

Prediction of the Secant Stiffness of the Yield Point
The same method is used to calculate the secant stiffness of the yield point. The formula is proposed in Equations (24)- (27).
A comparison of the test and parametric analysis results with the formula results is shown in Figures 49 and 50. The formula could fit the test with an error less than 20%. Most of the values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. The formula could accurately predict the secant stiffness of the yield point. A comparison of the test and parametric analysis results with the formula results is shown in Figures 49 and 50. The formula could fit the test with an error less than 20%. Most of the values of the FE model are higher than those of the formula, which indicates that the formula could conservatively predict the FE model results. The formula could accurately predict the secant stiffness of the yield point.

Comparison to Other Test Results
The formulas could accurately predict the ultimate strength capacity and stiffness of the test and parameter analysis accurately. To prove that the formula is universal for a similar composite shear wall, the formula is compared with the other results in the literature. Because the methods used for the stiffness of elastic and yield point are different, the elastic stiffness and secant stiffness of yield point are not compared. The comparison of the bearing capacities is shown in Table 2.

Comparison to Other Test Results
The formulas could accurately predict the ultimate strength capacity and stiffness of the test and parameter analysis accurately. To prove that the formula is universal for a similar composite shear wall, the formula is compared with the other results in the literature. Because the methods used for the stiffness of elastic and yield point are different, the elastic stiffness and secant stiffness of yield point are not compared. The comparison of the bearing capacities is shown in Table 2.

Comparison to Other Test Results
The formulas could accurately predict the ultimate strength capacity and stiffness of the test and parameter analysis accurately. To prove that the formula is universal for a similar composite shear wall, the formula is compared with the other results in the literature. Because the methods used for the stiffness of elastic and yield point are different, the elastic stiffness and secant stiffness of yield point are not compared. The comparison of the bearing capacities is shown in Table 2. The errors between the formulas and tests in other studies are not more than 20%, indicating that the formula in this paper is also suitable for predicting the ultimate strength capacity and yield strength capacity of composite shear wall in different crosssection forms. Thus, the calculated formula has a certain universality.

Conclusions
In this paper, FE models are developed to predict the hysteresis behaviour of CWSC. The mechanical mechanism of the wall is studied. The effects of key design parameters on the bearing capacity and stiffness are investigated. The main conclusions could be drawn as follows: 1. The web plate and concrete are the main components that resisted the lateral force.
The web plate is found to contribute to between 55% and 85% of the total shearing resistance of the wall. 2. The corner of wall mainly resisted the vertical force and the rest of wall resisted the shear force. The concrete is separated into several columns by the stiffened plates, each of which is independent and resisted by vertical force. 3. The elastic stiffness and ultimate strength capacity are enhanced with increasing wall thickness, steel ratio, axial compression ratio, and channel length-to-width ratio. The elastic stiffness and ultimate strength capacity are weakened with increasing shear 1.
The web plate and concrete are the main components that resisted the lateral force. The web plate is found to contribute to between 55% and 85% of the total shearing resistance of the wall.

2.
The corner of wall mainly resisted the vertical force and the rest of wall resisted the shear force. The concrete is separated into several columns by the stiffened plates, each of which is independent and resisted by vertical force. 3.
The elastic stiffness and ultimate strength capacity are enhanced with increasing wall thickness, steel ratio, axial compression ratio, and channel length-to-width ratio. The elastic stiffness and ultimate strength capacity are weakened with increasing shear span ratio. 4.
The parameters affecting the ductility of CWSC are the steel ratio, shear span ratio, axial compression ratio and the length-to-width ratio of the channel. The influence of steel ratio and shear span ratio on ductility is positively correlated, while the axial compression ratio and the length-to-width ratio of the channel are negatively correlated.

5.
Formulas are proposed to evaluate the ultimate strength capacity, yielding bearing capacity, elastic stiffness, and secant stiffness of the yield point of the composite shear wall. The formulas in this paper could predict the ultimate strength capacity more accurately than the formulas in specifications. Meanwhile, the formulas could predict the ultimate strength capacity in other tests from the literature well. The formulas can provide a basis for engineering design. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.